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^ ABSTRACT 

Star forming galaxies are thought to dominate the sub-mJy radio population, 
q but recent work has shown that low luminosity AGN can still make a significant 

5— 1 contribution to the faint radio source population. Spectral indices are an important 

tool for understanding the emission mechanism of the faint radio sources. We have 
d observed the extended Chandra Deep Field South at 5.5 GHz using a mosaic of 42 

pointings with the Australia Telescope Compact Array (ATCA). Our image reaches 
I an almost uniform sensitivity of ~12 /iJy rms over 0.25 deg 2 with a restoring beam 

j>. of 4.9 x 2.0 arcsec, making it one of the deepest 6cm surveys to date. We present the 

\Q 5.5 GHz catalogue and source counts from this field. We take advantage of the large 

C^) amounts of ancillary data in this field to study the 1.4 to 5.5 GHz spectral indices of 

CO the sub-mJy population. For the full 5.5 GHz selected sample we find a flat median 

^sO spectral index, a mcc j = —0.40, which is consistent with previous results. However, the 

qq spectral index appears to steepen at the faintest flux density levels (Ss.s&Hz < 0.1 

mJy), where a me d = —0.68. We performed stacking analysis of the faint 1.4 GHz 
£Nj selected sample (40 < Siaghz < 200 /iJy) and also find a steep average spectral 

7— I index, a = —0.8, consistent with synchrotron emission. We find a weak trend of 

steepening spectral index with redshift. Several young AGN candidates are identified 
• — h using spectral indices, suggesting Gigahertz Peaked Spectrum (GPS) sources are as 

common in the mJy population as they are at Jy levels. 

5_j 

C$ Key words: galaxies: evolution — radio continuum: galaxies 



1 INTRODUCTION 

One of the most fundamental issues in astrophysics is when 
and how stars, galaxies and black holes form, and how they 
evolved with cosmic time. The black hole mass and bulge re- 
lation ( Kormendy fc Richstone|l9 95 Magorrian et al.|1998| ) 
suggests there is a link between active galactic nuclei (AGN) 
and star formation, hence studying the interaction between 
stars, galaxies and AGN phenomena is crucial for under- 
standing galaxy formation in the early universe and how 
those galaxies evolve to the objects we see today. Radio 



E-mail: minh.huynh@uwa.edu.au 



emission can be produced by both AGN and star-forming 
processes, and thus radio wavelengths provide a unique win- 
dow to study the cosmic evolution of these two important 
processes. 

Bright radio sources (> 100 mJy) are associated with 
AGN activity (e.g. |Condo"n]|1984| [(Teorgakakis et al.||1999| |, 
but the Euclidean-normalized radio source counts flatten 
below about 1 mJy and this cannot be explained by a 
population of radio-loud AGNs. A non-evolving popula- 



tion of local (z < 0.1) low-luminosity radio galaxies Wall 
et al. ( 1986 1, strongly evolving normal spirals (|Condon|1984 



19891, and star forming galaxies ( jWindhorst et al.||1985 
Rowan-Robinson et al.|1993 1 have all been suggested to ex- 
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plain this new population. The faint radio source counts have 



been successfully modelled by star forming galaxies (Sey- 
mour et al.||2004 Huynh et al.|2005 l. The most commonly 
accepted paradigm has been that the sub-mjy population is 
largely made up of starforming galaxies. However a growing 
number of studies are finding that low luminosity AGN, both 
radio-loud and radio-quiet, still make a significant contribu- 
tion to the sub-mjy population Oarvis & Rawlings 2004 



2 THE OBSERVATIONS 
2.1 Observing Strategy 

We observed the extended Chandra Deep Field South at 6cm 
in two runs using the new Compact Array Broadband Back- 



Huynh et al.|2008||Seymour et al.|2008||Smolcic et al.|2008 
Padovani et al.|2009l 12011b 



There is conflicting evidence on the nature and prop- 
erties of faint radio sources, and in particular their spectral 
index properties and whether there is a flattening of the 



average spectral indices for faint radio sources (e.g. Ran- 
dall et al.p012| ). |Prandoni eFaT] pOOrl} found sources with 
1.4 GHz flux densities less than a few mjy had an aver- 
age 1.4 GHz to 5 GHz spectral index flatter than that of 



brighter radio sources. Owen et al. (2009) found a flattening 



of the average 325 MHz to 1.4 GHz spectral index below 
Siaghz < 10 mjy, but spectral indices appeared to steepen 
again at the faintest flux densities. In the Lockman Hole, 
the deepest radio field to date, no flattening of the 610 MHz 
to 1.4 GHz spectral indices was observed for Siaghz > 0.1 
mjy ( |Ibar et al.|2"009| >. 

Galaxies dominated by star formation processes are 
expected to have a spectral index of a = —0.8 (S oc 
u a , |Condon| |1992h , consistent with synchrotron emission 
from electrons accelerated by supernovae. Although ther- 
mal Bremsstrahlung (free-free) emission found in HII re- 
gions can have a flatter spectral index, a flat or inverted 
spectrum is usually attributed to the superposition of dif- 
ferent self-absorbed components of varying sizes at the base 
of the radio jet of a radio-loud AGN. A flattening in the 
average spectral index would imply there is a population 
of sources at sub-mjy flux densities with flat or inverted 
spectra, which are likely to be AGN. Studying the spectral 
index properties of the faint radio population is important 
for understanding the z — a relation, which is used to iden- 
tify the highest redshift radio sources (e.g.|De Breuck et al.| 



end (CABB; |Wilson et aL||20lT| ) with the full CABB 2048 
MHz bandwidth centred at 5.5 GHz. At 6cm the primary 
beam of the ATCA observations is ~10 arcmin (FWHM). A 
total of 42 pointings in a hexagonal layout (Figure [TJ was 
used to uniformly image the full 30 x 30 arcmin eCDFS 
area. Each pointing is separated by about 5 arcmin, which 
is approximately 0.5 FWHM of the primary beam at 6cm. 
The resulting mosaic is centred on the original GOODS field 
with approximate RA and Dec (J2000) of 3h32m20s and 
-27°48'34". 

The first run in August 2009 consisted of an initial set 
of 12 and 8 hour observations. A further 144 hours was al- 
located in January 2010. The ATCA was in the 6D and 6A 
configurations in August 2009 and January 2010, respec- 
tively, providing a maximum baseline of 6km. The mosaic 
pointings were cycled through rapidly with each pointing 
observed for only 1 minute in each pass. Taking into ac- 
count overheads such as calibrator scans, telescope drive- 
time and down time from inclement weather, the total ef- 
fective integration time of the mosaic is about 115 hours, 
or 2.7 hours per pointing. The secondary calibrator, 0347- 
279, was observed in the middle of the mosaic cycle for gain 
and phase calibration. During the August 2009 observations, 
secondary calibrator observations were separated by about 
20 minutes. The January 2010 observing run had uncharac- 
teristically poor phase stability and so the phase calibrator 
was observed about every 10 minutes instead. Primary flux 
density calibration was performed through observations of 
the source PKS B1934-638, which is the standard primary 
calibrator for ATCA observations (Ss.sghz ~ 5.1 JjQ. 



2.2 Data Reduction 

We used the Multichannel Image Reconstruction, Image 



2004||Klamer et al.|2006| ), identifying young radio AGN such Analysis and Display (MIRIAD) software package pault 



as Gigahertz Peaked Spectrum sources and Compact Steep 
Spectrum sources ( |0'Dea|[T998l |Randall et al.||2011j ), and 
determining whether star formation or AGN processes are 
responsible for the radio emission in these sources. 

Here we present 5.5 GHz observations by the Australia 



et al. 1995 I to reduce the CABB data. The latest version 



Telescope Large Area Survey (ATLAS, Norris et al. 2006 1 



team of the extended Chandra Deep Field South which are 
well-matched in area to the extensive mult i wavelength data 
in the region. This survey is important as it is the largest 
6cm survey deeper than 100 /xJy. We present the 5.5 GHz 
source catalogue and counts, and discuss the spectral index 
properties of the ATLAS radio sources. The paper is orga- 
nized as follows: In Section 2 we describe the observations 
and data reduction. We discuss the source extraction and 
present the source catalogue in Section 3. Source counts are 
derived from the catalogue and presented in Section 4. We 
perform a detailed spectral index (ctiACHz) analysis, and 
identify ultra-steep and Gigahertz Peaked Spectrum sources 
in Section 5. Concluding remarks are given in Section 6. We 
assume a Hubble constant of 71 km s" 1 Mpc -1 , J7m = 0.27 
and Qa = 0.73 throughout this paper. 



of MIRIAD was used to ensure tasks were updated to deal 
with data from this new back end. The MIRIAD task atlod 
was run with the "birdie" option to remove self-inteference 
from the 640 MHz clock harmonics and 100 channels at both 
band edges, where the bandpass response is poor. 

The 42 pointings were individually reduced and imaged. 
Automated flagging was performed with the task rmr/?a<j^] 
using an rms cutoff mode. In this mode RFI was identified 
if the rms in a channel exceeds the rms in a good channel 
by a set factor, or cutoff. This cutoff was 5a for the bright 
calibrators and 3cr for the field observations. To make sure 
the flagging was successful the data were inspected visually 
using the 3D visualisation tool plotvis 2 and any remaining 
RFI was manually excised. We explored a range of robust 



1 See ATCA Calibrator list: |http: //www. narrabri.atnf.csiro.au/ | 
calibrators/ 

2 http: //www. atnf.csiro.au/people/Emil.Lenc/tools/Tools/ 
Mirflag_Plotvis.html 
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Right Ascension (J2O00) 



Figure 1. LEFT: Greyscale of the full eCDFS 5.5 GHz mosaic. The 42 mosaic pointing centers are marked by crosses. The inner red 
rectangle marks the area covered by the MUSYC opt ical imaging ||Taylor et al,||2009| and the outer black rectangle marks the area 
covered by Very Large Array (VLA) 1.4 GHz imaging | |MiIler et al.|2008[ l. RIGHT: A cutout of a central portion of the mosaic to show 
more detail. 



weights in the imaging stage (Briggs 1995) as well as uni- 
form and natural weighting. We found a robust weight of 1 
resulted in a good compromise between sensitivity and res- 
olution. An image 2300 x 2300 pixels in size with 0.5 arcsec 
pixels was generated for each pointing to ensure recovery of 
the full field of the primary beam. Cleaning was performed 
using multi-frequency cleaning task mfclean with the clean 
region set to 8 arcmin from the center, which extends be- 
yond the 10% response level of the primary beam and hence 
encompasses almost all the area of interest. One iteration 
of self-calibration was performed for each pointing and the 
images were restored with a beam of 4.9 x 2.0 arcsec, the 
average synthesized beam of the 42 pointings. The MIRIAD 
task linmos was then used to stitch together the individual 
pointings to create a single mosaic map, shown in Figure [l] 
The CABB bandwidth is so large that we also split the 
2000 channels to make four sub-band images of approxi- 
mately 500 channels each. The central frequencies of these 
sub-band images are 4.80, 5.25, 5.72 and 6.19 GHz. We use 
these images to examine the SED of our sources and identify 
Gigahertz Peaked Spectrum sources in Section 5.3. A high 
resolution image with a beam of 2.9 x 1.0 arcsec was also 
generated with super-uniform weighting, and this is used to 
examine radio morphology ( Dehghan et al.|2011 l. 



2.3 Image Analysis 

The primary beam correction at the high and low end of the 
2 GHz band differs by about 50% at the half-power point of 
the primary beam for 5.5 GHz. This may affect the recovered 
flux density for sources at the edge of a pointing since linmos 
uses a single primary beam correction. To check this effect 
we ran simulations for sources with spectral indices varying 
from a = 1 to a = —2 (S oc v a ). We find the flux density 
of a flat spectrum source is recovered to within 1% in the 2 
GHz band at the half-power of the primary beam. The flux 



density is recovered to ~ 3% for a canonical synchrotron 
spectrum of a = —0.8. It is only for the steepest spectrum 
sources (a = —2) that the flux density is over-estimated, 
by ~ 7%, at the half-power point of the primary beam, and 
we don't expect these sources to be common. We therefore 
conclude that as long as the source lies within the half-power 
point (i.e. within 0.5 FWHM of a pointing center) this effect 
is minimal. 

We examined the noise characteristics of the mosaic 



using the Sextractor software package (Bertin & Arnouts 



19961. While SExtractor was developed for the analysis of 



optical images, it is known to work well on radio images also 



I Bondi et al.|2003|[Huynh et al.|2005l|Prandoni et al.|2006 1 

The choice of mesh-size is important. A mesh-size which 
is too small would result in the estimation being affected 
by individual sources, while larger mesh sizes can miss sys- 
tematic small-scale variations in the noise. Previous surveys 
have found that mesh-sizes with widths of 8 to 12 times the 
synthesised beam, produce good results for noise estimation 
in deep radio continuum surveys ( Huynh et al.|2005 Schin 
ncrcr et al.|p007|[2010l ). A more detailed analysis of simu- 
lated radio continuum images confirms that a mesh-size of 
10 to 20 synthesized beams produces optimum background 



noise estimates (Huynh et al. 20121. To estimate the local 



background noise variations we therefore used a mesh-size 
of ~ 10 x 10 synthesized beams (67 x 67 pixels). The SEx- 
tractor generated noise image is shown in Figure [2] (left), 
limited to areas within 5 arcmin, ~0.5 FWHM of the pri- 
mary beam. We find that the mosaic has an rms noise of 
~11.9 /iJy/beam, but with a large tail at higher noise levels 
(Figure [5] right). This tail is the result of (i) higher noise re- 
gions around extended sources that aren't cleaned perfectly 
and (ii) the greater noise at the edge of the mosaic, which 
are furthest away from any one pointing and therefore have 
the greatest primary beam correction. 

The effective integration time of the mosaic is approxi- 
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Figure 2. LEFT: Greyscale of the noise image from SExtractor, generated using a mesh-size of 10 X 10 synthesized beams. The inner 
red rectangle marks the area covered by the MUSYC optical imaging ( |Taylor et al.|2009[ l and the outer black rectangle marks the area 
covered by VLA 1.4 GHz imaging (Miller et al.i.2008}. RIGHT: The pixel distribution of the SExtractor noise image. Solid black line 
indicates the distribution for full noise image on the left, and the red dotted line indicates the distribution in the MUSYC imaging area. 
The peak of the distribution is at ~ 12 /^Jy. The high noise tail is mostly at the edges of the mosaic, where the primary beam response 
is lower. 



mately 115 hours. According to the online ATCA sensitivity 
calculator, ATCA should reach 6 /iJy/beam rms in 12 hours 
at 5.5 GHz with the new CABB upgrade. This corresponds 
to a sensitivity of 6 x \/42/l-4 = 28 /xJy/beam rms for a 
nyquist-sampled 42 pointing mosaic in 12 hours. We there- 
fore should have achieved . 30 ~ 9.7 /^Jy/beam. The 

a/115/12 

sensitivity we achieve is within ~ 20% of this value esti- 
mated from the ATCA sensitivity calculator. 



3 SOURCE EXTRACTION 

There are many radio source identification and measurement 
tools in common use, including MIRIAD/AIPS Gaussian 



which are parameters found by Huynh et al. (2012) to re- 



fitting routines imsad, sad and vsad, sfind (Hopkins et al 



2002[), Ducham p ( |Whiting||2012[ ) and SExtractor ( |Bertin fc 
Arnouts|1996 1 . Most of these use a simple S /N thresholding 
technique whereby a source with a peak flux density, or pixel 
value, above a set threshold (usually some multiple of the lo- 
cal noise) is deemed to be a true source. The false-discovery 
rate (FDR) approach is a robust statistical procedure which 
compares the distribution of image pixels to that of an image 
of equal size containing only noise, and the user threshold 
then places a limit on the fraction of identified sources which 
may be false, based simply on the statistics of the noise dis- 
tribution (IMiller et al.||2001|j. The FDR method has been 



implemented for source measurement in radio images in the 
MIRIAD task sfind ( |Hopkins et al.|2002[ ). The FDR method 
was found to be more reliable than simple S/N threshold 
methods for parameters which return similar completeness 
( Huynh et al.|2012[ >, so we therefore use sfind for the source 
detection. 

We restricted the search area to within 5 arcmin (~ 0.5 
FWHM of the primary beam) of the outer mosaic point- 
ing to minimise primary beam effects, and ran sfind with 
'rmsbox' set to 10 synthesized beams and 'alpha' set to 1, 



turn reasonable noise estimates, completeness and reliabil- 
ity. Setting the 'alpha' parameter to 1 would return a list of 
sources which is 99% reliable for the case of a perfect image 
with pure Gaussian noise, sfind returned an initial list of 
143 sources for inspection. Each source was then individu- 
ally characterised as both a point source and Gaussian using 
the MIRIAD task imfit. 

Ten radio sources show classical core-lobe AGN mor- 
phology, and hence are clearly components of a single source 
(see Figure |3|. These sources were fitted as multiple Gaus- 
sians with imfit and the multiple components are listed indi- 
vidually in the final catalogue. In total there are 123 sources 
and 142 source components in the final source list. 



3.1 Deconvolution 

The ratio of the integrated flux to the peak flux is a direct 
measure of the extension of a radio source: 



"Stot/ -Speak — 



mm^max 



/ ^min^n 



(1) 



where # m i n and # m i n are the source FWHM axes and 6 m i n 
and 6max are the synthesized beam FWHM axes. This ratio 
can be used to determine whether a source is resolved. 



An analysis similar to Prandoni et al. ( 2006 1 and Huynh 



et al. ( 2005 ) was performed to determine which sources are 



resolved. Whether a source is successfully deconvolved de- 
pends on the S /N ratio of the source and not just the image 
beam-size. Using the Gaussian fits, we examined the ratio of 
integrated flux density to peak flux density as a function of 
source signal to noise (Figure Eb. Assuming the sources with 
•Stot /-Speak < 1 are due to noise then the upper envelope is 
defined as 



Stot/Speak = 1 + 10/(1 + [Speak/a]) 1 



(2) 
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Figure 4. The ratio of integrated (Stot) fl ux density to peak flux 
density (Speak) a- s a function of source signal to noise (Speak/ 17 ) • 
The dotted line shows the upper and lower envelopes of the flux 
ratio distribution that contains 90% of the unresolved sources. 
The large dots indicate sources which are deconvolved successfully 
and considered resolved. 



This envelope is found by determining a lower curve that 
contains 90% of the Stot /Speak < 1 sources (see Figure [4| 
and then mirroring it on the Stot /Speak > 1 side. Only 
47/143 (33%) sources lie above the upper envelope and are 
considered to be successfully deconvolved (i.e. resolved). 



3.2 The Source Catalogue 

The source catalogue is reported in Table [I] Point-source 
measurements are given for sources which are not success- 
fully deconvolved. The fitted peak, integrated source flux 
density and deconvolved source sizes from the Gaussian 
fits are given for the resolved, or successfully deconvolved, 
sources. Absolute calibration errors dominate for high S/N 
sources, but internal fitting errors shown in Table [l] domi- 
nate for the majority of sources, which are low S/N. 

Column (1) - ID. A letter, such as 'a', 'b', etc., indicates 
a component of a multiple source. 

Column (2) - Source IAU name 

Columns (3) and (4) - Source position: Right Ascension 
and Declination (J2000) 



Column (5) - Peak source flux density (/iJy) 

Column (6) - Uncertainty in peak source flux density 

(MJy) 

Column (7) - Integrated flux density (/^Jy). Zero indi- 
cates source is not successfully deconvolved and hence no 
integrated flux density is given. 

Column (8) - Resolved flag. Zero indicates source is not 
successfully deconvolved. 

Column (9) - Deconvolved major axis (arcsec). Zero in- 
dicates source is not successfully deconvolved. 

Column (10) - Deconvolved minor axis (arcsec). Zero 
indicates source is not successfully deconvolved. 

Column (11) - Deconvolved position angle (degrees), 
measured from north through east. Zero indicates source is 
not successfully deconvolved. 

Column(12) - Local noise level, rms, in /xjy. 

3.3 Comparison with VLA survey 

The CDFS area has been observed with the Very Large Ar- 
ray (VLA) at 4.9 GHz. Four VLA pointings were used to 
cover a region of approximately 20 x 20 arcmin in size, with 
depths varying from 7/iJy/beam rms at the pointing centers 
to 50 /xJy/beam rms at the edges ( Kellermann et al.|2008 |. 
The resolution of the VLA image is approximately 3.5 arc- 
sec, which is almost the same size as the synthesized beam 
of our ATCA imaging. 

We compared the single component sources in our cat- 
alogue with those in the VLA catalogue. We find 49 of our 
sources have a VLA 4.9 GHz flux density measurement from 
(2008). A straight comparison of the flux 



IKellermann et al.l 



densities (Figure [5 left) suggests that ATCA flux densities 
are higher than the VLA flux densities, and the ATCA/VLA 
flux density ratio has a mean of 1.28 ± 0.09 and median 
of 1.14. However one source with Ss.sghz ~ 12.2 mjy is 
a clear outlier. Subsequent analysis finds this source is an 
AGN with QSO-like optical colours and hence variability 
could be an issue with this source. Excluding this outlier, 
the ATCA/VLA flux density ratio has a mean of 1.21 ± 
0.06 and median of 1.13. For a spectral index of a = —0.8 
we expect the ATCA flux densities to be about 8% less than 
the VLA measurements, if the VLA and ATCA are cali- 
brated on the same scale. To minimise spectral index effects 
we compare the VLA flux densities to those of our sources 
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Table 1. The ATLAS 5.5 GHz catalogue 



ID IAU name R.A Dec Speafc dS pea fc Si n ± Dccon Dccon Dccon Dccon & local 

(J2000) (J2000) (>Jy) (M-TyJ <>Jy) Flag Bmajor Bminor PA 



3 


ATCDFS5 


J033348 


.75-280232 


.8 


3 


33 


:48 


.75 


-28:02 


:32 


.8 


354 


.8 


69. 


6 





.0 








.00 





.00 





.0 


55 


6 


4 


ATCDFS5 


J033346 


59-273405 


.4 


3: 


33 


:46 


.59 


-27:34 


:05 


.4 


181 


.5 


37. 








.0 








.00 





.00 





.0 


46 


9 


5 


ATCDFS5 


J033341 


32-273809 


.5 


3 


33 


:41 


.32 


-27:38 


:09 


.5 


216 


.9 


22. 


2 





.0 








.00 





.00 





.0 


36 


.8 


6 


ATCDFS5 


J033338 


.35-280031 


.0 


3: 


;33 


:38. 


.35 


-28:00 


: 3 1 





536 





20. 








.0 








.00 





.00 





.0 


27 


.2 


7 


ATCDFS5 


J033334 


59-274751 


.2 


3 


;33 


:34 


.59 


-27:47 


: 5 1 


.2 


148 


.1 


16. 


5 





.0 








.00 





.00 





.0 


23 


6 


8 


ATCDFS5 


J033333 


43-275332 


.8 


3 


;33 


:33 


.43 


-27:53 


:32 


.8 


467 


.7 


14. 


1 





.0 








.00 





.00 





.0 


18 


.4 


9 


ATCDFS5 


J033332 


57-273539 


.2 


3 


;33 


:32 


.57 


-27:35 


:39 


.2 


460 


.6 


21. 


8 





.0 








.00 





.00 





.0 


19 


.5 


10 


ATCDFS5 


J033327 


53-275726 


.6 


3: 


;33 


:27 


.53 


-27:57 


:26 


.6 


109 


.4 


15. 


1 





.0 








.00 





.00 





.0 


16 


.5 


11 


ATCDFS5 


J033325 


.85-274343 


.2 


3 


33 


:25 


.85 


-27:43 


:43 


.2 


186 


.2 


19. 


9 





.0 








.00 





00 





.0 


15 


.0 


12 


ATCDFS5 


J033322 


75-275500 


.1 


3 


33 


:22 


.75 


-27:55 


00 


.1 


100 


.9 


17. 


1 





.0 








.00 





.00 





.0 


15 


1 


13 


ATCDFS5 


J033321 


31-274138. 


.6 


3 


33 


:21 


.31 


-27:41 


:38 


.6 


267 


.8 


18. 


4 





.0 








.00 





.00 





.0 


13 


.4 


14 


ATCDFS5 


J033318 


71-274940 


.0 


3 


;33 


:18 


.71 


-27:49 


:40 


.0 


72 


. 1 


12. 


2 





.0 








.00 





.00 





.0 


13 


.7 


15 


ATCDFS5 


J033318 


30-273440 


.4 


3 


;33 


: 18 


.30 


-27:34 


:4() 


.4 


96 


.3 


12. 


9 





.0 








.00 





00 





.0 


14 


.2 


16 


ATCDFS5 


J033316 


96-274121 


.7 


3: 


;33 


:16 


.96 


-27:41 


:21 


.7 


78 


.3 


14. 


1 


125 


.2 


1 


2 


.82 


1 


12 


-50 


.1 


11 


.7 


17 


ATCDFS5 


J033316 


76-280016 


.2 


3 


33 


:16 


.76 


-28:00 


:16 


.2 
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The ATLAS 5.5 GHz survey of the Extended Chandra Deep Field South 7 



Table 1. continued 
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as measured in the ATCA 4.8 GHz sub-band image (Fig- 
ure^ right). We find that the ATCA/VLA ratio, excluding 
the outlying source, has a mean of 1.23 ± 0.08 and median 
of 1.18. The ATCA flux densities are therefore about 20% 
greater than the ones measured by the VLA. 

Discrepancies in measured flux densities can result from 
the details of the fitting algorithm used in the source extrac- 
tion, but this is expected to affect only faint sources. The 
higher ATCA flux density is clear in the strong sources with 
0.3 < Ss.sghz < 1 mjy, not just in the faint sources, which 



suggests that choice of fitting algorithm is not the cause of 
the differences. A possible discrepancy between the VLA and 
ATCA flux density scales has been noted in the past ( |Norris| 



et al. 2006 Kellermann et al. 20081, but with ATCA flux 



densities reported to be lower than VLA ones by 5 to 20%. 
However, those comparisons were done at 1.4 GHz, and so 
it is possible that the calibrator scales assumed can lead to 
different flux density ratio comparison results between 1.4 
GHz and 5 GHz. 
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Figure 5. LEFT: Comparison of the flux densities in Table[T]with corresponding measurements in Kelle rmann et ah] ( |2008| ) . The dotted 
line denotes ATCA flux density = VLA flux density. The ATCA values appear greater (by ~20%) than the VLA ones, but there is 
considerable scatter at the lowest flux densities. RIGHT: Same as the left, but the comparison is with ATCA flux densities measured in 
the 4.8 GHz sub-band image. 



3.4 Completeness and Flux Boosting 

The completeness of the source catalog was estimated by 
simulations. Monte-Carlo simulations were performed by in- 
jecting 8000 sources at random locations of the map and ex- 
tracting them using the same technique as adopted for the 
production of the catalog. While completeness levels vary 
significantly across the image, becoming worse at the edges 
of the mosaic due to primary beam effects, we can recover 
the overall completeness level of the generated catalog by 
injecting sources over the full image. A single source was 
injected per simulation, and the input flux density varied 
from 30 to 3000 /iJy to sample the full range of interest. 
The completeness as a function of flux density is shown in 
Figure [6] (left). The completeness rises steeply from about 
10% at 50 /iJy to approximately 90% at 130 ^Jy. The 50% 
completeness level occurs at approximately 75 /^Jy. 

Source detection algorithms that rely on finding a peak 
above a local noise background can lead to flux boosting. 
This is because sources that lie on a noise peak have their 
flux density increased and therefore have a higher proba- 
bility of being detected, while sources which lie on a noise 
trough have decreased flux densities and may be excluded 
altogether. This effect is most pronounced in the faintest 
flux density bins, i.e. in sources with a SNR close to the 
limit of the catalog. The degree of flux boosting can be esti- 
mated by examining the output to input flux density of the 
simulations (Figure [6j right). We find that flux densities are 
boosted by about 5% to 1% for measured flux densities of 
100 to 200 /iJy, on average. The flux boosting is negligible 
for sources with true flux densities brighter than 200 fiJy. 

The simulations can also be used to gauge the positional 
accuracy of the catalog by comparing input and output posi- 
tions. The median of the RA and Dec offsets are plotted for 
various input flux density bins in Figure [7] The positional 
accuracy can be estimated from the standard deviation in 
the offsets of each of these bins. We find that at the faintest 



levels (50 fJ.3y) the fitted RA and Dec uncertainties are ap- 
proximately 0.24 and 0.46 arcsec, respectively. The total po- 
sitional accuracy is ~0.3 arcsec or better for sources that are 
0.1 mjy and brighter. 

3.5 Source Size and Resolution Bias 

Weak extended sources with large total integrated flux den- 
sities may have peak flux densities that fall below the detec- 
tion threshold. To derive source counts which are complete 
in terms of total flux density the so-called resolution bias 
must be determined. We follow the formalism of P randonil 



et al. ( 2001 1 and Huynh et al. ( 2005 I in calculating the res- 



olution bias. 

The maximum size (8 mal ) a source of total flux density 
Stot can have is determined by 



Stot/Cde 



'max 



/ ^min^n 



(3) 



where b m i n and & max are the synthesized beam FWHM axes 
and Odet is the detection limit. As the SFIND detection limit 
varies across the image, we use the 50% completeness level, 
as determined by the simulations of Section |3.4| The re- 
sulting # max as a function of total flux density is plotted in 
Figure [8] (dotted line). We plot the angular sizes (0) of the 
catalogued sources as a function of total flux density in Fig- 
ure [8] also. The angular size 6 is defined as the geometric 
mean of the fitted Gaussian major and minor axes. We find 
that the largest catalogued sources are in good agreement 
with the # max function. 

The minimum angular size (#min) is estimated from 
Equation 2, with a equal to the noise of the full image (12 
^iJy) . The resulting # m i n as a function of total flux density is 
plotted in Figure[8] (solid line). This 6 m ia constraint is impor- 
tant at low flux density levels, where # max becomes unphys- 
ical (smaller than a point source). Overplotted in Figure [8] 
(dashed lines) is the expected median angular size obtained 
from Windhorst et al. 1990 relations for a 1.4 GHz sample: 
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Figure 6. LEFT: Completeness as a function of input flux density, as derived from the Monte-Carlo simulations. Completeness is 
the number of extracted sources divided by number of input sources. RIGHT: The distribution of (input flux density - output flux 
density) /input flux density as a function of output flux density for the simulated sources. The solid red line is the median of the 
simulation and the dashed lines mark the 1 sigma upper and lower bounds. The effect of flux boosting at the faint end is dramatically 
illustrated by the rapid downturn to negative values below about 0.1 mjy. 
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Figure 7. The offset in RA (left) and Dec (right) between the recovered positions of sources in the simulation and the true input 
positions, as a function of input flux density. The error bars mark the 1 sigma uncertainty in the position as a function of input flux 
density. 



a _ o" c°-30 

fmed — ^ >J 1.4GHz 

where 5*i.4ghz is in mjy. The extrapolation to 5.5 GHz was 
done assuming a spectral index of 0, -0.5 and -0.8 between 
1.4 and 5.5 GHz. The integral angular size distribution, h(8), 
is | | Windhorst et al.|1990[ ): 



h(6) = exp" 



-ln2(9/e mod )° 



At the bright end our source sizes are consistent with the 
|Windhorst et al.] ( |1990[ ) relation, however most of the sources 
are unresolved and therefore we can not draw any conclu- 
sions about the full sample. 



As outlined in Prandoni et al. (20011, the overall an- 
gular size limit, #u m = max(# max , # m i n ) and the expected 
integral size distribution, h(9), allows an estimation of the 
fraction of sources larger than the maximum detectable size, 
and hence missed by the survey. The resolution bias correc- 
tion factor is then simply ■ This correction factor is 
plotted in Figure [9] (left). It has a maximum of about 1.35 
at a flux density of 100 /^Jy, where the limiting overall angu- 
lar size, #ii m , becomes dominated by (9 m in- Windhorst ct al. 



rection to be about 10%. The caveat to the resolution bias 
correction is that the integral angular size distribution at 
these low flux densities is not well known. High resolution 
MERLIN and VLA observations of 92 radio sources with 
Siaghz > 40/iJy ( |Muxlow et al.||2005| ) hint that the faint 
radio population may have an integral angular size distribu- 
tion different to the brighter sample. The integral angular 



size dist ribution for the|Muxlow et al. ( |2005p sample differs 
from the |Windhorst et al.| ( |1990[ ) relation by up to 20%(Fig- 
ure[9j right). High resolution observations are likely to miss 



low surface brightness galaxies, but Muxlow et al. (20051 



estimate that they have missed only 10% of the faint radio 
sources in their field. If the Muxlow et al. ( 2005 1 sizes are 



( 1990 1 estimate the uncertainty in the resolution bias cor- 



used the resolution correction factor near 100 /iJy would be 
approximately 5%, and zero elsewhere. 



4 SOURCE COUNTS 

The differential radio source counts were constructed from 
the 6cm eCDFS catalogue of Section 3. In computing the 
source counts we used integrated flux densities for extended 
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Table 2. The 5.5 GHz 


source 


counts. 
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(MJy) 






(sr Jy ) 


(X10 ) 
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1.68 X 10 10 


0.47 ± 0.14 
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Figure 8. The fitted angular size as a function of total flux den- 
sity. Sources fitted or deconvolved as a point source are plotted 
with an angular size of zero. The solid line indicates the minimum 
angular size (0 m j n ) of sources in the survey, below which decon- 
volution is not considered meaningful. The dotted line shows the 
maximum angular size (0 max ) above which the survey becomes 
incomplete due to resolution bias. The dashed lines indicate the 
median source sizes expected from the Windhorst et al. 1990 re- 
lation, as a function of flux density, for a spectral index of 0, -0.5 
and -0.8 between 1.4 and 5.5 GHz. 



sources and the point source fit for all other sources. The 
components of multiple sources were summed and counted 
as a single source. For comparison with with other 6cm 
studies, the source counts are normalised to a non-evolving 
Euclidean model. At 6cm the standard Euclidean integral 
counts are N(> Se cul ) = 60 x Sq c ]^ sr _1 , where Seem is in 



Jy ( |Fomalont et al.||l99l| Ciliegi et al.||2003| >. The results 
are summarised in Table |2J where for each bin we report 
the flux density interval, median flux density, the number of 
sources detected (N), the number of sources after complete- 
ness and resolution bias corrections have been applied (Nc), 
the differential source count (dNc/dS), and the normalised 
counts (Nc /-/Vexp) • The counts are normalised to -/V exp , the 
number expected in the bin from the standard Euclidean 
count. The Poissonian errors in the count are 
where C is the total correction factor, Nc/N. The estimated 
total error in the counts is the Poissonian error with the res- 
olution bias uncertainty (10%) and completeness correction 
uncertainties (2 - 4 %) added in quadrature. 

Our results are compared with previous work in Fig- 
ure |10| In this Figure we only show source counts below 
4 mjy. The counts are well known above this flux density 



(Altschuler 1986 Wrobel & Krause 1990 Prandoni et al 



20061: there is an initial steep rise between 10 and Uy, an 



excess with respect to the Euclidean prediction between 1 
Jy and 0.1 Jy and a steep slope in the counts for the flux 
density range 100 to 0.5 mjy. Some authors claim there is a 



flattening of the 6cm source counts at mjy levels (e.g. Ciliegi 
et al.|[2003 l, while others do not see the flattening down to 
0.5 mjy flPrandoni et al.|2006 | 



Our source counts are consistent with other published 
6cm source counts for Seem > 0.3 mjy. At fainter flux den- 
sities our source counts are significantly lower than |Ciliegi| 
et al. ( 2003 1 counts by a factor of about two. Our counts are 
consistent with Donnelly et al. ( |1987[ ) and Fomalont et al. 
1991, however the counts in the faintest bins (S < 50/iJy) 



from Fomalont et al. 1991 are about a factor of two higher 
than our counts in the lowest bin. Fomalont et al. 1991 cata- 
logued sources to about 4a in their image, so it is likely that 
they have spurious sources in their faintest bins. We note 
that out survey area is almost 3 times greater than Ciliegi 
et al. 2003 (0.25 deg 2 versus 0.087 deg 2 ) and 5 times greater 
than F omalont et al. 1991 (0 . 25 de g 2 versus 0.05 deg 2 ). Us- 
ing the Driver & Robotham ( 2010 1 expression for quantify- 



ing cosmic variance, we find that our survey has a cosmic 
variance of about 35 to 90% for bins of dz — 0.1 in size at 
z < 1. Similarly the Ciliegi et al. 2003 survey has a cosmic 
variance of 49 to 115%. The difference in the counts at the 
faint end can therefore be attributed to cosmic variance. 

Finally we note the 6cm surveys in the literature have a 
central frequency of 5 GHz and the difference of 0.5 GHz in 
the observing frequency will have an impact on the source 
counts. To estimate the counts at 5 GHz we applied spectral 
indices of 0, -0.4 and 0.8 to our catalogue and recalculated 
the source counts. We find that the source counts change 
by about 5% - 40% but this does not account for the sig- 
nificantly lower counts seen in the faintest bins of our data 
compared to previous 6cm surveys. 



5 SPECTRAL INDEX ANALYSIS 

5.1 1.4 to 5.5 GHz Spectral Indices 

To investigate the spectral index properties of faint radio 
population we matched the 5.5 GHz catalogue to sources 
in the Very Large Array (VLA) 1.4 GHz survey of the ex- 
tended Chandra Deep Field South ( |Miller et al.|2008| |. This 
survey is a good match in the area covered, as seen in Figure 
1. The depth of the dMiller et al.||2008l) survey is about 8 



/iJy rms, which should be sensitive enough to detect coun- 
terparts to all but the faintest 5.5 GHz sources with flat or 
inverted spectra. Importantly, the beam of the VLA obser- 
vations is 2.8 X 1.6 " beam, which is only a factor of ~1.5 
better than our observations. With similar resolutions these 
images have a similar surface brightness sensitivity, and thus 
the measured flux densities can be used directly for spectral 
index analyses. 

To begin the spectral index analysis we removed the 
multi-component sources from the analysis as their inter- 
pretation is complicated by the core-jet structure, leaving 
113 individual 5.5 GHz sources for investigation. 101/113 
sources have a 1.4 GHz match within 2 arcsec (FWHM of 
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Figure 9. LEFT: The resolution bias correction as a function of flux density, assuming the Windhorst et al. 1990 integral angular source 
size distribution. RIGHT: The integral angular size distribution from the Windhorst et al. 1990 relation (for flux densities as shown) and 
that observed in the Muxlow et al. 2005 sample. The Muxlow et al. 2005 sample has a median flux density of 75 ^tJy. 
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Figure 10. Normalized 5.5 GHz differential source counts for different samples: Prandoni et al. (2006) (empty blue squares); Ciliegi et 
al. (2003) (empty circles); Fomalont et al. (1991) (red upside-down triangles); Donnelly et al. (1987) (empty triangles). The eCDFS 5.5 
GHz source counts presented in this work (filled squares) are corrected for completeness and resolution bias as explained in the text. 
Vertical bars represent Poisson errors on the normalized counts. 



the synthesised beam of the VLA observations). The rest 
were inspected and we found four were out of the 1.4 GHz 
image area, and seven had a faint counterpart in the 1.4 
GHz image that was below the Miller et al. (2008) cata- 
logued threshold. The 1.4 flux density for these sources was 
measured manually with the MIRIAD task imfit. In sum- 
mary 108/109 (99%) of the 5.5 GHz sources in the 1.4 GHz 



image area have a 1.4 GHz counterpart, and hence a spectral 
index measurement. The source without a 1.4 GHz counter- 
part, ID82, has a 7a limit of Siaghz < 48pJy, and hence 
is inverted with a > 0.55. The median spectral index for 
the ATCA 5.5 GHz sample is a mcd = —0.40 (see Figure 



111, which is consistent with that of the deepest 6cm sam- 
ple to date (Q m cd = —0.38 for Ss.sghz > 16/iJy, Fomalont 
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Figure 11. Spectral index distribution for sources in the ATCA 
6cm sample. The vertical dashed line indicates the median value 
of the sample (a me£ j = —0.40). 



ct al. 1991). Figure 12 (left) shows the spectral index as a 



function of 5.5 GHz flux density. The median spectral in- 
dex is a mc d = —0.35 for 0.1 < Ss.bghz < 1 nijy. This is 
consistent with published values for 6cm selected sources of 
similar flux density, for example Prandoni et al. 2006 who 
found Q mc d = —0.4 for 0.4 < S^ghz < 4 mjy and Donnelly 
et al. 1987 who found a mcd = -0.42 for 0.4 < S 5G Hz < 1.2 
mjy. 

It has been claimed that spectral indices flatten from 
bright (mjy) to fainter (sub-mjy) flux densities. For ex- 
ample, Ciliegi et al. 2003 observed a flattening of the ra- 
dio spectral indices in 6cm sources with a mo d = —0.37 for 
0.1 < Ssghz < 0.2 mjy and a mcd = -0.81 for S 5G Hz > 0.2 
mjy. Our median spectral index for 0.1 < Ss.sghz < 0.2 
mjy is —0.37, which is consistent with this result, however 
we do not have enough mjy sources to confirm the spec- 
tral flattening from the mjy to sub-mjy regime. |Owen et al.| 
( |2009[ > also observed the spectral flattening from mjy to sub- 
mjy flux densities, albeit between 325 MHz and 1.4 GHz, 
but they found that the radio spectra steepened again at the 
faintest flux densities (Siaghz < 0.1 mjy). At the faintest 
flux density levels in our sample, Ss.sghz < 0.1 mjy, the 
median spectral index is a mc d = —0.68 showing a trend to 
steeper spectra at the faintest flux densities. This is a direct 
confirmation of what |Owen et al.| ( |2009[ ) concluded through 
stacking analysis. 

The fraction of flat or inverted sources (a > —0.5) is 
55% for the whole sample, consistent with the findings of 
a high proportion (50 - 60%) of flat or inverted sources in 
/xjy samples selected at similar frequencies (Fomalont et al 



T991) |Windhorst et al.|p93] |Prandoni et aL|2006[ ), but the 

fraction of flat or inverted sources drops to 11/30 (36%) for 
Ss.sghz < 0.1 mjy. The synchrotron-like median spectral 
index for Ss.sghz < 0.1 mjy is consistent with star-forming 
galaxies dominating the population at these faintest flux 
density levels. 

We also examined the spectral indices of the 1.4 GHz se- 
lected sample. Of the 514 1.4 GHz sources in the Miller ct al. 
(20081 observations, 3 are not in our 5.5 GHz catalogued re- 



gion. We found 36 1.4 GHz sources that were associated with 
a multiple 6cm source and these were also removed from 
further spectral index analysis. 106/475 (22%) of 1.4 GHz 



sources have a 6cm counterpart within 2 arcsec (Figure 12 
right). The sfind a parameter we used is best approximated 
by 5a for the noise characteristics of our image. Upper limits 
to the spectral index assuming a 5<r detection limit at 5.5 



GHz are thus plotted in Figure 12 (right) to illustrate where 



the spectral index information becomes incomplete. Figure 



12 (right) shows that the 1.4 GHz sources start to be missed 
in the 5.5 GHz image for Siaghz < 0.5 mjy. 

To derive median spectral indices of the full 1.4 GHz 
sample at faint flux densities we performed a stacking analy- 
sis in bins of 1.4 GHz flux density. The 6cm data was stacked 
at the positions of 1.4 GHz sources not detected in the 6cm 
image. The effective pixel value in the stack was taken to 
be the median value of the stacked positions, and the flux 
density determined by a point source fit to the signal at 
the centre of the stack. Offset stacks were generated by ran- 
domly offsetting each stacked position by up to 8 arcsec. 
The uncertainty in the flux density of the median stack is 
then estimated as the standard deviation of the flux density 
measured in 300 random offset stacks. The uncertainty from 
the point source fit is then added in quadrature to this es- 
timate to derive a total uncertainty in the median stacked 
flux density. The results are summarised in Table [3] As a 
check of our stacking technique, we also derived the noise- 
weighted flux densities from the stacks and found agreement 
(to within 9%) with the median stacked flux density. 

The median spectral index for 1 < Siaghz < 10 mjy 
is a — —0.57, consistent with the Prandoni et al. (2006) 1.4 
GHz sample. The stacking results show that the spectral 
index of the faint (40 < Siaghz < 200 /iJy) sources range 
from a = -0.64 ± 0.22 to -0.92 ± 0.08, consistent with 
synchrotron emission from star forming galaxies. Again, this 
is consistent with Owen et al. 2009 who find radio sources 
are steep at the faintest flux density levels (Siaghz < 100 
pjy). We find that 17/22 (77%) of faint (Siaghz < 200 
/xjy) sources detected at both frequencies are flat or inverted 
(a > —0.5), but stacking analysis of the full sample of faint 
/xjy level 1.4 GHz sources shows they are steep on average 
and flat or inverted sources do not dominate at the faint flux 
density levels. 

Finally we note that the differences between VLA and 
ATCA flux densities (Section 3.3) implies a possible system- 
atic effect on the measured spectral indices. If the ATCA 5.5 
GHz flux densities are actually over-estimated then the true 
spectral indices are steeper by 0.13. 

5.2 Optical counterparts and redshifts 

Optical counterparts for our 5.5 GHz sample were identi- 
fied from the COMBO-17 (Wolf et al. 200 4) and Multi- 
wavel ength Survey by Yale-Chile (MUSYC, |Taylor et "aL 
120091) imaging surveys of the eCDFS. Both surveys are about 
30 x 30' in size and hence well matched in area to our 5.5 
GHz survey. Spectroscopic redshifts are available for many 



sources in the eCDFS (e.g. Vanzella et al. 2008 Balestra 



et al.|[2010 Mao et al. in prep) and we searched these cat- 



alogues for redshifts of the radio sources, adopting photo- 
metric redshifts only when spectroscopic redshifts are not 
available. We prioritize MUSYC photometric redshifts over 
COMBO-17, due to the larger filter set used. As in Sec- 
tion 5.1, we exclude the multiple sources from this analysis. 
We find 104/113 sources lie within MUSYC or COMBO-17 
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Figure 12. LEFT: 1.4-5.5 GHz spectral index vs. 5.5 GHz flux density for the ATCA 5.5 GHz sample. Only sources catalogued at 5.5 
GHz are shown. The dashed line indicates the median spectral index, a = —0.40. RIGHT: 1.4-5.5 GHz spectral index vs. 1.4 GHz flux 
density for the |Miller et al.| | |2008| l 1.4 GHz sample. Only sources catalogued at 1.4 GHz are shown (black dots), and black arrows mark 
upper limits assuming a detection limit of 5cr in the 5.5 GHz image. The average spectral indices from stacking analysis of the faint 1.4 
GHz sources are shown as red squares. The dashed red line indicates the median spectral index of the 1 < SiagHz < 10 m Jy sample, 
a = —0.57. The black dotted line denotes the canonical spectral index from synchrotron emission, a = —0.8. 



Table 3. Results of stacking 1.4 GHz positions in the 5.5 GHz image. N^t is the number of 1.4 GHz sources in each bin which are 
detected in the 5.5 GHz image. N s t a ck is the number of 1.4 GHz sources not detected in the radio image and therefore the number 
stacked in this analysis, resulting in the stacked S^^GHz listed in Column 5. Stacked a is the spectral index for the stacked sample. 
Average a is the spectral index for all sources in the flux density bin, obtained by combining the detections with the stacked results. 
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(MJy) 




N det 
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mean SiagHz 

(MJy) 


stacked S 5 . 5G Hz 

(MJy) 


stacked a 


average a 


40 < SiAGHz 


< 60 


1 


19 


52.0 ± 2.8 


21.6 ± 6.4 


-0.64 ± 0.22 


-0.59 ±0.22 


60 < SlAGHz 


< 80 





83 


71.4 ± 1.6 


25.8 ± 3.7 


-0.74 ± 0.11 


-0.74 ± 0.11 


80 < SiAGHz 


< 100 


2 


102 


89.0 ± 1.5 


29.5 ± 2.2 


-0.80 ±0.06 


-0.79 ±0.06 


100 < SlAGH 


» < 130 


4 


76 


114.4 ± 1.9 


32.4 ± 3.7 


-0.92 ±0.08 


-0.88 ±0.08 


130 < SiaGH 


, < 200 


15 


66 


156.7 ± 2.1 


48.3 ± 5.2 


-0.86 ±0.08 


-0.75 ±0.09 



imaging, of which 97/104 (92%) have counterparts within 
2 arcsec (~ 0.67 FWHM of the synthesized beam). Spec- 
troscopic redshifts are available for 45 sources. We adopt 
MUSYC photometric redshifts for 47 sources and COMBO- 
17 photometric redshifts for 5. 

The radio luminosities and redshift distribution for 
these sources are shown in Figure |13| The median redshift of 
the 5.5 GHz sources is z — 0.87, and the distribution has a 
tail to high redshifts. The survey is sensitive to sources with 
luminosities as faint as Ls.sghz = 10 23 W Hz -1 up to z ~ 1, 
indicative of low powered AGN or powerful starbursts. At 
higher redshifts the sources have AGN-like radio powers. 



5.3 The z-a Relation and Ultra-Steep Spectrum 
Sources 

Here we further explore the spectral index properties of our 
sample, with a focus on the extent to which the observed 
correlation of radio spectral index with redshift, the z — a 
relation, persists for the faint flux-density population. The 
z— a relation is one of the most successful tracers of high red- 
shift radio galaxies (HzRGS) flTelens et al.||1979| |Blundell 
et al.||1999| |Klamer et aT1|2006| |Ker etal.||2012^ Although 
an ultra-steep spectrum (USS, a < — 1) doesn't guarantee 
a high redshift radio source, a higher fraction of high red- 



shift sources can be found amongst USS sources (Roettger- 
ing et al.|19"97l|De Breuck et al.|2000||2001| . Until recently, 
USS samples have been limited to Siaghz > 10 mjy since 
tens or hundreds of square degrees are required to identify 
a meaningful sample. A fainter (SewUBz > 100 ^3y) pop- 
ulation of 58 USS sources was found in the Lockman Hole. 
Optical/IR counterparts for ~60% of this sample suggests 
they have redshifts spanning 0.1 < z < 2.8, and the non- 
detection of the other ~40% suggest they lie at even higher 
redshift ( jAfonso et al.|2011 i. 



Figure [14| (left) shows a versus redshift. While a linear 
fit shows only a weak correlation (a correlation coefficient of 
r = —0.17, which only becomes weaker if the highest red- 
shift data point is excluded, r = —0.13), there is, however, a 
clear excess of steeper spectra at z > 2. While there are very 
few systems at z > 2, reflecting the sparseness of HzRGs on 
the sky and the difficulty in identifying them with small area 
surveys such as this, the fact that 5/8 of the sources at z > 2 
have a < —0.7 is suggestive that the z — a relation persists 
even to the faint flux-densities measurable here. To ensure 
that this result is not simply a consequence of Malmquist 
(right) shows that there is essentially no 



14 



bias, Figure 

correlation between luminosity and spectral index (a cor- 
relation coefficient of r = —0.08), demonstrating that we 
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Figure 13. LEFT: The 5.5 GHz radio luminosity as a function of redshift. These are k-corrected assuming the spectral indices measured 
between 1.4 and 5.5 GHz. RIGHT: The redshift distribution of the radio sources. 



are not biased against flat or inverted spectrum systems at 
high-z. 

Our results are in agreement with those for much 
brighter samples ( |Ker et al.|2012|. Also c onsistent with the 
Lockman Hole results from |Afonso et aL (2011 1, our steep- 
est spectrum sources have redshifts spanning 0.5 < z < 2.7. 
We do, however, have a much higher completeness of op- 
tical counterparts for this population (~ 80% compared to 
~ 60%), even when including a further six 1.4 GHz sources 
not detected at 5.5 GHz that have a < —1.0. This suggests 



that, consistent with Afonso et al. ( 2011 1, the USS technique 



for selecting high redshift radio galaxies remains efficient to 
sub-mjy levels, although the small areas currently accessible 
with deep surveys will limit the numbers of such rare high 
redshift objects able to be identified. 

5.4 Candidate Young AGN Identified Through 
Spectral Indices 

One approach to identifying candidate young AGNs is 
through their size and spectral index properties. Radio 
sources with compact sizes and steep spectra (Compact 
Steep Spectrum, CSS) or peaked spectra (Gigahertz Peaked 
Spectrum, GPS), are potential candidates for an early stage 
of the AGN evolutionary process ( |Fanti et al.|1995||Tinti~fe 
de Zotti 20061. At bright radio flux-densities such sources 



have sizes comparable with or smaller than the optical host 



galaxy (O'Dea 19981, but still show double-lobed struc- 
tures on small scales (~ 1 — lOkpc) (O'Dea 1998 Snellen 
|et al.|1998| . An alternative hypothesis to explain their small 
sizes is that, rather than being young, they are 'frustrated' 



sources, being confined by a dense medium (O'Dea et al. 
|1991[ ). CSS and GPS sources are relatively common in bright 
source samples, accounting for up to 30% of sources selected 



at 2.7 and 5.0 GHz (e.g. |KapahipMl||Peacock fc Wall|1982 1. 
A recent unbiased sample selected at 2.7 GHz finds more 
than 10% of bright (Jy) sources are CSS or GPS sources 
(Randall et al. 20111. Existing CSS and GPS samples are 



largely heterogeneous, however, and there are few samples 



at faint flux densities (mjy versus Jy level, e.g. Randall et al 



20121 



In total, we identify five GPS candidates that are con- 
sistent with being either a young AGN or an AGN confined 



by a dense medium, based on their size and spectral shape 
(full details are given in Appendix |Aj. More complete radio 
spectral information is required to obtain a complete sam- 
ple of GPS sources in this field, but our preliminary analysis 
shows at least 5/40 sources with S5.5GHZ > 0.5 mjy meet 
the GPS criteria. This GPS fraction is similar to that found 
in Jy level samples ( jRandall et al.||2~011| ), suggesting GPS 
sources remain an important population to mjy levels. 



6 CONCLUDING REMARKS 

We have presented new observations at 5.5 GHz of the ex- 
tended Chandra Deep Field South. The 0.25 deg 2 region was 
observed with the Australia Telescope Compact Array us- 
ing a mosaic of 42 pointings. The resultant image reaches 
an almost uniform noise level of ~12 fiJy rms and has a 
resolving beam of 4.9 x 2.0 arcsec. Using a false-discovery- 
rate method, we extracted 123 sources or 142 source com- 
ponents. Ten sources were resolved multiple sources with 
jet-lobe structure and hence fitted as multiple components. 

After carefully correcting for completeness, flux boost- 
ing and resolution bias, we derived source counts at 5.5 GHz. 
These are amongst the deepest source counts at 6cm but 
from an area 3 to 5 times larger than the previous counts 
to these depths. The ATLAS 5.5 GHz counts are consistent 
with the counts derived from other 5 GHz surveys at brighter 
flux densities, but are lower than counts in the literature by 
a factor of two for S^.sghz < 0.3 mjy. This discrepancy is 
attributed to cosmic variance because of the small size of 
the surveys involved. This fluctuation in the 5.5 GHz source 
counts at the faint end is similar to that seen at 1.4 GHz for 
SiagHz < 0.1 mjy ( jNorris et al.|2011 l 

The 1.4 - 5.5 GHz spectral index has also been deter- 
mined for all the sources in this field. We find a median spec- 
tral index for the ATCA 5.5 GHz sample of a mc d = —0.40. 
We find the median spectral index for the faintest flux den- 
sity levels in our sample, Sz,.zghz < 0.1 mjy, is a mc d = 
—0.68. The 1.4 GHz sources start to be missed in the 5.5 
GHz image for Siaghz < 0.5 mjy so we performed a stack- 
ing analysis to examine the average spectral index of the full 
faint 1.4 GHz selected sample. From the stacking analysis we 
find that the faintest 1.4 GHz sources (40 < Siaghz < 200 
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Figure 14. LEFT: The spectral index, a, as a function of redshift. There is a weak trend of decreasing a with increasing redshift. The 
dotted line shows the linear best fit. RIGHT: The spectral index, a, as a function of 5.5 GHz radio luminosity. The dotted line shows 
the linear best fit, but there is no trend of a with radio luminosity. 



/xJy) also have steep spectra. This is consistent with the re- 



sults of Owen et al. ( 2009 1 who found that radio spectra 
steepens at the faintest flux density levels (Siaghz < 100 
pijy), after flattening from mjy to sub-mjy flux densities. 

The 5.5 GHz sources have been cross-matched to exist- 
ing optical/NIR data to obtain redshifts and we find, con- 
sistent with earlier work, that the z — a relation seems to 
persist to these low flux densities. In addition, several candi- 
date young AGN have been identified using the peak in the 
spectral energy distribution between 1.4 GHz and 6 GHz, 
suggesting that GPS sources are as common in the mjy 
population as they are at Jy levels. 

In future work we will explore the spectral index prop- 
erties of the ATLAS sources against other attributes, such 
as polarization and X-ray hardness, with the aim of distin- 
guishing low-luminosity AGN from starforming galaxies in 
the sub-mjy population and determining how the observed 
properties of the sub-mjy population are linked to these two 
physical emission processes. 
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APPENDIX A: THE GPS CANDIDATES 

We combined the four sub-band images from our 6cm obser- 
vations with the VLA 1.4 GHz observations to identify GPS 
source candidates. Based on the spectral index properties 
alone, we identify 3 strong GPS candidates, and another 4 
weaker candidates where there is possibly a peak but the 
large errors in the flux densities makes the existence of the 
peak less clear. The SEDs for these sources are shown in 
Figure [AT] and |A2"| ID64 is included in the weak category as 
the 1.4 GHz flux density is greater than the 6cm sub-band 
GHz flux densities, and there is a possibility that the SED 
may not peak at all in between 1.4 and 6.2 GHz. 

Subsequent measures of source size show that all but 
two sources have compact (< 10 kpc) radio morphology. Of 
the two extended sources, ID131 is a classical FRII double 
with lobes separated by about 21 kpc and ID64 is a FRI 
radio galaxy with lobes extending out to about 110 kpc. 
These two sources are rejected as GPS candidates on the 
basis of their extent, and in the case of ID64 also on the 
spectral shape. 

Optical counterparts are taken from the MUSYC and 
COMBO-17 data, and the properties of the GPS candidates 
are summarised in Table [ATI The radio luminosities are in 
the observed frame, i.e. not fc-corrected, as the spectral index 
between 1.4 and 5.5 GHz is uncertain for these sources. 

In summary, all our low flux-density GPS candidates 
are consistent with FRI or FRII radio sources, at redshifts 
about 2 ~ 0.5 and higher. They are dominated by those with 
physical extents around 5-10 kpc, well within the character- 
istic size of massive galaxies at these redshifts. As such, they 
are likely to be young jets yet to have pushed outside the 
stellar disks of their host galaxies. This is consistent with the 
properties of GPS sources identified at the brightest radio 
flux densities. 

We discuss the individual GPS candidates in detail be- 
low. 

ID85 

GPS candidate ID85 is a strong (~ 12 mjy) point 
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Figure Al. The radio spectral energy distributions of the 3 strong GPS candidates. The red square marks the total S§,$ghz from the 
full CABB band image and the other 6cm points are from the four CABB sub-band images. The 1.4 GHz point is from Miller et al. 2008. 
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Figure A2. The radio spectral energy distributions of the 4 weak GPS candidates. The red square marks the total Ss.sGHz from the 
full CABB band image and the other 6cm points are from the four CABB sub-band images. The 1.4 GHz point is from Miller et al. 2008. 



source at 5.5 GHz. It is only marginally resolved in the 
higher resolution 1.4 GHz data, with a deconvolved size of 
0.7 x 0.5 arscec. The bright (Rvcga = 18.63) optical coun- 
terpart appears to be a point source and has QSO colours. 
The COMBO-17 photometric redshift of z = 1.57 ± 0.09 is 
consistent with the spectroscopic redshift of z — 1.5542. At 
this redshift the deconvolved angular size implies this radio 
source is about 6 kpc in extent and has a radio luminosity 
of Lxaghz = 2.48 x 10 25 W Hz -1 . Hence this source has 
a radio power consistent with FRI radio galaxies, or low- 
powered FRIIs. 
ID93 

This S5.5GHZ = 3.4 mjy source is not resolved at 5.5 
GHz and only marginally resolved at 1.4 GHz. The bright 
(Rab = 18.84) optical counterpart appears to be a point 
source and has QSO colours. The MUSYC photometric red- 
shift of z = 0.55 ± 0.01 is consistent with the spectroscopic 
redshift of z — 0.5444. At this redshift the deconvolved an- 
gular size implies this radio source is about 10 kpc in extent. 
The radio luminosity of Liagh* = 1.25 x 10 24 W Hz -1 is 
consistent with a low luminosity FRI or a very powerful 
star burst. Higher resolution radio imaging is required to es- 
tablish whether the extended emission shows a morphology 
characteristic of jets (most likely) or of diffuse star forma- 
tion (less likely) . If the radio emission was all associated with 
star formation, it would imply a star formation rate (SFR) 
of ~ 100 M Q yr~\ which is high for galaxies at this redshift, 
but not extreme. 

ID 131 

This bright (Ss.sghz = 12.10 mjy) source is extended 
at 5.5 GHz with a deconvolved size of 3.5 x 1.0 arcsec. An 



inspection of the 1.4 GHz image shows a compact double, so 
we manually refit the 1.4 GHz source as a double Gaussian, 
finding a total 1.4 GHz flux density of 39.5 mjy. The optical 
counterpart in MUSYC is a Rah ~ 22.01 galaxy with a 
photometric redshift of z — 1.75 ± 0.01. The lobes in the 
1.4 GHz image are ~2.5 arcsec apart, which corresponds to 
a linear extent of 21 kpc. The radio luminosity of L\aghz 
= 3.06 x 10 26 W Hz -1 is consistent with this source being a 
FRII, or high-powered FRI, radio galaxy. 
ID64 

This bright (S$.$ghz = 12.10 mjy) source is extended 
with a deconvolved size of 3.4 x 0.9 arcsec at 5.5 GHz. It 
is also extended at 1.4 GHz and shows an FRI-like mor- 
phology with faint lobes. This source is not easily fit as a 
multiple Gaussian. A manual fit to the 1.4 GHz core of the 
source yields a 15 mjy source with deconvolved size of 3.2 x 
1.3 arcsec. The bright (i?v cg a = 19.00) optical counterpart 
appears to be a point source, has QSO colours, and a possi- 
ble nearby neighbour 5 arcsec to the east. The COMBO-17 
photometric redshift of z — 0.75 ± 0.15 is consistent with 
the spectroscopic redshift of z — 0.7386. The optical source 
to the east does not have a reliable photometric redshift in 
either catalogue, possibly due to another nearby faint ob- 
ject to its east which could be contaminating the source's 
photometry, so it is unclear whether a merger is seen in the 
optical bands. At z = 0.7386 this radio source has a linear 
extent of about 23 kpc, but the faint FRI lobes are about 15 
arcsec in distance or ~ llOkpc apart. The radio luminosity 
of Liaghz = 2.21 x 10 25 W Hz -1 is consistent with that of 
an FRI radio galaxy. 

ID 92 
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Table Al. Summary of the properties of the GPS candidates. The 3 strong candidates are listed before the weak candidates. The radio 
luminosities have not been k-correctcd. 



ID 


SlAGHz 


3$. 5GHz 


R mag 


phot z 


spec z 


angular-size 


linear-size 


LlAGHz 


L5.5GHZ 




(mjy) 


(mjy) 








(arcsec) 


(kpc) 


(W Hz" 1 ) 


(W Hz" 1 ) 


85 


4.00 


12.25 


18.63 a 


1.57 ± 0.09 a 


1.5542 6 


0.7 X 0.5 


6.0 x 4.3 


2.48 x 10 25 


7.61 x 10 25 


93 


1.65 


3.37 


18.84 c 


0.55 ± 0.01 c 


0.5444 6 


1.2 x 0.9 


10.0 x 5.7 


1.25 x 10 24 


2.56 x 10 24 


131 


39.5 d 


12.10 


22.01 c 


1.75 ± 0.01 c 




2.5 e 


21.4 e 


3.06 x 10 26 


9.38 x 10 25 


64 


15.3 f 


6.74 


19.00 a 


0.75 ± 0.15 a 


0.7423 i > 


3.2 x 1.3 


23.3 x 9.5 


2.21 x 10 25 


9.72 x 10 24 


92 


0.265 


0.616 


21.63 c 


0.73 ± 0.01 c 


0.733 6 


1.2 x 0.6 


8.8 x 4.4 


3.72 x 10 23 


8.66 x 10 23 


100 


1.91 


2.39 


22.30 c 


0.81 ± 0.01 c 




0.8 x 0.3 


6.0 x 2.3 


3.29 x 10 24 


4.12 x 10 24 


103 


0.66 


0.79 


23.63 c 


1.57 ± 0.01 c 




0.8 x 0.4 


7.1 x 3.8 


4.18 x 10 24 


5.00 x 10 24 



a COMBO-17, Wolf et al. 2004, Vega-magnitude 
b ATLAS, Mao et al. 2012 

c MUSYC, Cardamone et al. 2010, AB-magnitude 

revised using a double Gaussian fit to Miller et al. (2008) 1.4 GHz image 
e distance between lobes in Miller et al. (2008) 1.4 GHz image 
* revised using a free Gaussian fit to Miller et al. (2008) 1.4 GHz image 

This source, with a 5.5 GHz flux density of 0.6 mjy, is 
the faintest GPS candidate. It is a point source in the 5.5 
GHz image but marginally resolved in the 1.4 GHz image 
with a deconvolved size of 1.2 x 0.6 arcsec. The source has 
an Rab = 21.63 counterpart with a photometric redshift 
of 0.73 ± 0.01, consistent with the spectroscopic redshift 
of z = 0.733. At this redshift the deconvolved angular size 
implies this radio source is about 9 kpc in extent and has 
a radio luminosity of L 1AG Hz = 3.72 x 10 23 W Hz" 1 . This 
radio power is indicative of a low-luminosity AGN or high 
star formation (SFR ~ 40 M yr -1 ). 

ID 100 

This S5.5GHZ = 2.4 mjy source is not resolved at 
5.5 GHz and only marginally resolved at 1.4 GHz. The 
Rab = 22.30 optical counterpart has a photometric redshift 
of 2 = 0.81 ± 0.01 but there is no spectroscopic redshift. At 
this photometric redshift the radio source is about 6 kpc in 
extent. The radio luminosity of Liaghz = 3.29 x 10 24 W 
Hz -1 is consistent with an FRI radio galaxy. 

ID 103 

This source is relatively faint for a GPS candidate 
(S5.5GHZ = 0.79 mjy). It is not resolved at 5.5 GHz and only 
marginally resolved at 1.4 GHz. The faint Rab = 23.63 opti- 
cal counterpart has a photometric redshift of z = 1.57±0.01 
but there is no spectroscopic redshift. At this photometric 
redshift the radio source is about 7 kpc in extent. The radio 
luminosity of Liaghz = 4.18 x 10 24 W Hz -1 is consistent 
with an FRI radio galaxy. 
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